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Abstract 

Three dimensional Thirring model with N four- component Dirac fermions, reformulated 
as a lattice gauge theory, is studied by computer simulation. According to an 8 3 data and 
preliminary 16 3 data, chiral symmetry is found to be spontaneously broken for N = 2, 4 and 
6. N = 2 data exhibits long tail of the non- vanishing chiral condensate into weak coupling 
region, and N = 6 case shows phase separation between the strong coupling region and the 
weak coupling region. Although the comparison between 8 3 data and 16 3 data shows large 
finite volume effects, an existence of the critical fermion flavor number N cr (2 < N CI < 6) for 
which the chiral behavior changes its character, is suggested by the current numerical data. 

PACS number(s): 11.15.Ha, 11.30.Qc, 11.30.Rd 



Thirring model, appeared in the (1 + 1) dimensional world around forty years ago, has 
shed light on understanding various intriguing issues of field theories, e.g., exactly-soluble 
models, bosonization, conformal fixed point, etc. Recently, the very model with N four- 
component Dirac fermions has attracted attention in the (2+1) dimensional world with 
regard to renormalizability of four-fermion models in the context of 1/iV-expansion fiJ-0, 
dynamical mass generation for the fermion in the analysis of Schwinger-Dyson equation 
0-0 and so on. In the listed works, usually 1/N expansion technique is employed since 
interesting features of the model are inherently non-perturbative. In this note, using a Monte 
Carlo simulation technique, we consider the Thirring model on a three-dimensional (3D) cubic 
lattice, and ask whether the vector-type four fermion interaction of the model breaks chiral 
symmetry spontaneously. 

The Lagrangian for the 3D Thirring model in Euclidean space time is given by 

G 

^ = E^«(7m^ + + ^(E^7^«) 2 - (!) 

a * a 

Typical way of treating this model is to introduce an auxiliary vector field ^ ^ f|, ^ 0] . 
On a lattice, this leads to an auxiliary vector field formulation [flO | . An alternative approach 
which is attractive to us, is to regard as a £7(1) gauge field. Indeed a reformulation of 



the Thirring model as a gauge theory can be achieved through a hidden local symmetry [11 



£hls = £&(7A + m)^ a + - VNd^) 2 , (2) 

where ip a is a 4-component Dirac spinor, is a fictitious Nambu-Goldstone boson field (or 
equivalently a Stiickelberg field), and = <9 M — -^A^. 

Noticing that the unitary gauge-fixed form of the above U(l) gauge theory coincides ex- 
actly with the original Thirring model, we expect benefits from gauge-invariant formulation. 
One such benefit is that, in the limit of zero bare fermion mass, one can avoid complications 
like the dynamical generation of parity- violating mass [|12| or Chern-Simons term ||13|| . The 
other is that we have a freedom to choose the most appropriate gauge for each particular 
situations of interest, e.g., analysis of the Schwinger-Dyson equation in a rigorous fashion || 
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and a construction of a gauge-invariant effective potential for the chiral order parameter |J. 

Since the model we are considering is free from parity violation by an exact argument 
in the gauge theory formulation, we can focus on the spontaneous breakdown of chiral 
symmetry (actually a flavor symmetry breaking of U(2N) — > U(N) x U(N)). There are 
various interesting aspects to be addressed once the above chiral symmetry breaking occurs, 
which were raised by analytic methods. First, a rather systematic approach of Schwinger- 
Dyson equations in the gauge theory formulation || suggests a finite, critical number of the 
four-component fermion species, 

N cr = ^| « 4.32 (3) 

which is obtained in the infinite limit of four-fermion coupling g (g is a dimensionless coupling, 
G rescaled by the ultraviolet cutoff A, g = GA). In contrast, the auxiliary field formulation 
with the "Landau gauge fixing 11 implies otherwise.^ A crude constant mass approximation 
to this formulation seems to predict a chiral symmetry breaking for all N || 

/7T 2 \ 

g cr oc exp(^— NJ, (4) 

which also tells us the ratio of the values of the critical coupling g CT for different iV's. Until 
now, there is no known equivalent relationship between g cr and A" in the gauge theory 
formulation. However, if such a relation exists, it should be different from the above equation 
since there exists a critical number of fermion flavor A" cr in the gauge theory formulation 
||14|| . Second, it may be interesting if we can read the chiral condensate (ipip) as a function of 
the coupling g for a given N. In the gauge theory formulation ||, a critical scaling behavior 
was suggested to follow 

(V>V) « a n exp ( — 7T=l) • 

Of course, on a finite lattice, the behavior of the chiral condensate will be modified by finite 
volume effect and may be difficult to compare with this formula away from g CT . 



1 In fact this Landau gauge fixing amounts to an gauge in the hidden local symmetry formulation if 
we concede that the auxiliary field is a gauge field. 
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Analytic studies in the framework of the Schwinger- Dyson equations combined with an 
1/N expansion raised the questions in the above and it is interesting to investigate these 
questions by a lattice study. In choosing an appropriate action for numerical simulation, 
one can use an auxiliary vector field formulation without gauge symmetry or employ a 
gauge field formulation based on a hidden local symmetry. Recently simulation results of 



the auxiliary field formulation case was reported in Ref.|T0[. Here we adopt the gauge 
theory formulation and take advantage of various benefits coming from the gauge theory 
formulation. A discretized version of the 3D action, Cl with a staggered fermion formulation 
for the fermion part, is 

£l = E $ t a(MtM)- 1 $ a -N/3J2 cos (<f>(x +fi)- <P(x) + O^x)), (6) 

a fi 

where $ is the pseudo fermion and 6^ is the gauge field, 

M x , y = \ E ^(*)(e"* (!e) <W - eT^-^Sy^) + m5 x , y (7) 

with staggered fermion phase factor, rj^x), and j3 = 1/g, after dropping the constant piece. 
We include a finite bare fermion mass term and investigate the massless limit by simulating 
successively smaller fermion mass, m, since our numerical method prohibits simulating the 
chiral limit directly. The continuum limit of this lattice model is free from parity-violating 



mass or Chern-Simons term [15|| . The gauge theory formulation of 3D Thirring model in 
Eq.(^|) looks like 3-dimensional quantum electrodynamics (QED) + scalar QED with kinetic 
term of the gauge field and the Higgs degrees of freedom truncated. Note that the link 
variable exp(i0 M ) can not be used for the vector-fermion interaction term, tpA^^ip, in the 
auxiliary field formulation. This was not encountered in lattice QED3 case because the 
gauge symmetry fixes uniquely the interaction between the fermion and the gauge field and 
it is the same as Eq.([7|). Since this theory is non-asymptotically free, such differences may 
lead to different continuum limits. In addition, since in our gauge field formulation the 
gauge kinetic term will be induced by quantum effects, we can not say whether it will follow 
compact or non-compact form yet. Such an effect seems to be realized indeed by some kind 
of discontinuous behavior in the actual simulations for the strong coupling region [|T0|. 



Thus far, we have completed simulations for N = 2, 4 and 6 systems on a 8 3 lattice 
volume with fermion mass, m = 0.05, 0.025, and 0.0125 by use of hybrid Monte Carlo 
algorithm (HMC). Preliminary data on 16 3 lattice for N = 2 and 6 will be used for finite 
volume effect study. As usual, by defining fermion fields on even sites only, one can avoid 
an additional fermion doubling from M^M and can simulate even N system with this HMC 
algorithm |17]]. For odd N's, we need to use other numerical algorithm such as hybrid 
molecular dynamics method flUfl . Typically, we have accumulated r ~ 10, 000 time unit after 
discarding initial r = 100 time unit as the equilibration time with choosing appropriate Ar(= 
0.025,0.02,0.0125) which gives reasonable acceptance rate (50% to 90%). The Metropolis 
step interval is set to r = 0.5. As an order parameter for the chiral symmetry, we calculate 
the chiral condensate 

(W) = (TrM" 1 ), (8) 

where we use a single noise vector to estimate the condensate. In addition, we also calculate 
internal energy (Ep) defined by (N(3(3 — ^ cos (<f>(x + /i) — <f>(x) +9 fl (x)))), and specific 
heat (Ep 2 ) — (Ep) 2 for each N. Mass spectra for the fermion, 7r(pseudo scalar channel) and 
p(vector channel) are also calculated. However, due to a short extent in the time direction, 
we were not able to fit the propagators reliably. Thus, in the following discussion, we will 
rely mainly on the chiral condensate measurement. 

Qualitative nature of our numerical simulation result can be summarized in Fig. 1. The 
figure is the m — > extrapolated chiral condensate on a 8 3 lattice volume for N = 2, 4, 
and 6 respectively as a function of (3, where we assumed a naive linear behavior in m for 
the chiral condensate (~ c + c\m). The figure is similar to QED3 case; in [fTTf] , N = 2 
shows smooth change of the chiral condensate and a long tail into weak coupling regime, 
N = 3 shows a complex behavior of the condensate near the critical point, and N = 5 shows 
an abrupt change of the condensate. This supports the existence of the critical number 
of fermion flavors for chiral symmetry breaking in QED3 |Tj|. Likewise, if we believe that 
the naive linear extrapolation of the condensate is a correct procedure in the gauge theory 
formulation of 3D Thirring model, there appears to be a critical N around which the phase 
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transition changes its character. A complicated second order phase transition which connects 
the strong coupling region to the weak coupling region smoothly (N = 2) turns into a phase 
transition which has a phase separation (N = 6). 
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Figure 1: The chiral condensate, {ipip) m ^ , vs. f3 for iV = 2 (O points), iV = 4 (+ 
points), and N — 6 (□ points) on a 8 3 lattice. 
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Figure 2: j3 2 (tpip) m ^ vs. /3 for A" = 2 on a 8 3 lattice. 

To elucidate these behaviors more clearly, in Fig. 2, we search for a continuum scaling 
window for N = 2. Like QED3 N = and 1 case, the chiral condensate exhibits a long 
plateau in weak coupling region. On the other hand, N = 4 shows two slopes and hints flat 
region. 
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Figure 3: (V^)m->o vs. (3 for iV = 2(0), = 4(+), and = 6(D) on a 8 3 lattice. 



In Fig. 3, we plot (V ; V ; )m-»o vs - m or der to see whether the critical behavior follows a 
mean field theory, 

(W)~{p c -0)h. (9) 

It appears that N = 2 data does not agree with a mean field behavior and = 6 data 
seems to follow a mean field expectation. For A^ = 4 data, it is very hard to say whether 
the system behaves in one way or the other. From internal energy density and specific heat 
peak, we obtain (3 C = 0.9(1) for N = 2, 0.35(5) for N = 4, and 0.18(1) for N = 6. These 
values are roughly consistent with those obtained from Fig. 1, 2, and 3. 

According to the QED3 case, large finite volume effect is expected ]T7|, [B|. In order to 
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purchase this effect, simulations on larger lattices and smaller m values are needed. Thus, 
we are investigating the same system on a 16 3 lattice with the masses m = 0.05, 0.025 and 
0.0125. In Table 1, we compare {ipip) m _^ for N = 2 and in Table 2 compare that for 
N = 6 on 8 4 and 16 3 volume. The first column is the chiral condensate for 8 3 using a 
linear extrapolation of all range of m data, the second column is that for 16 3 using a linear 
extrapolation of smaller m data (m = 0.025, 0.0125), and the third column is that for 16 3 
using a quadratic fit to existing data. Although the dependence of the chiral condensate 
on both the size of lattices and the extrapolation method are considerable quantitatively 
(particularly in the very weak coupling region of N = 2 case), we found that the overall 
qualitative nature of the chiral condensate does not seem to be changed: In all three cases, 
N = Q data are compatible with mean-field behavior but N = 2 data is not. 



p 


<^>(8 3 ) 


(^)(16 3 ) 


<^>«w(16 3 ) 


0.5 


0.589(3) 


0.641(9) 


0.643 


0.75 


0.460(2) 


0.493(8) 


0.488 


1.0 


0.142(2) 


0.207(4) 


0.195 


1.25 


0.0218(9) 


0.053(2) 


0.048 


1.3 


0.0134(8) 


0.019(2) 


0.007 


1.4 


0.0071(7) 


0.003(2) 


-0.009 



Table 1: Comparison of the chiral condensate for N = 2 
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(^>(8 3 ) 


(^)(16 3 ) 


(^) a w(16 3 ) 


0.18 


0.135(2) 


0.263(4) 


0.239 


0.2 


0.017(1) 


0.085(2) 


0.058 
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-0.004(1) 


0.031(2) 


0.017 


0.25 


0.0014(8) 


0.004(9) 


-0.003 
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0.0027(7) 


-0.028(1) 


-0.010 


0.35 
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0.003 


0.4 


0.0004(6) 


0.0004(8) 


-0.003 


0.45 




0.0006(10) 


-0.001 


0.5 


0.0002(5) 


0.0005(7) 


-0.001 



Table 2: Comparison of the chiral condensate for N = 6 



In conclusion, we have found spontaneous chiral symmetry breaking for N = 2, 4 and 6 



in 3D Thirring model by use of Monte Carlo simulation on a 8 3 and a 16 3 lattice. However, 
the chiral condensate in the m — > limit suggests that the nature of the phase transition in 
N = 2 is different from that in N = 6. N = 2 case has a long tail of non- vanishing chiral 
condensate in weak coupling region and N = 6 case has a phase separation between the strong 
coupling region and the weak coupling region. It suggests a possibility that there exists N CT 
where the chiral phase transition becomes a mean field type transition. This situation is 
quite similar to QED3. In contrast to the lattice simulation result of the auxiliary vector 
field formulation, our result seems to be more consistent with the gauge theory formulation 
of 3D Thirring model in continuum. 

For a future work, we think that further detailed study on finite volume effects is needed 
although we have discussed our preliminary data on a 16 3 lattice. Also, two further questions 
remain open: One is whether lattice 3D Thirring model is similar to 3D lattice QED. The 
other is whether the induced gauge field kinetic term takes the form of compact QED-like or 
noncompact QED-like. The latter issue will be determined by the importance of topological 
sector. At this moment, we do not know what the induced term looks like. 
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